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A non trivial application of a modern computer language (" C" ) in a highly structured and 
object-oriented fashion is presented. The contest is that of Lie algebra representations 
(irreps), specifically the problem of reducing the products of irreps with the weight tree 
algorithm. The new WBase 2 . version with table-generation and Young tableaux display 
capabilities is introduced. 

1. Introduction 

Calculations in algebra representation theory, in particular decompositions of prod- 
ucts of irreps, are needed in several sectors of physics. The Dynkin approach to the 
representation theory^ is known by physicists thanks to its generality: all simple Lie 
algebras, included the exceptional ones, are described and manipulated in the same 
formal environment. In the Dynkin approach the algebras are described uniquely 
by the I x I Cartan matrix, where I is the rank of the algebra. The irreps of a given 
algebra are identified by a unique highest weight vector of / positive integers. 

The purpouse of this contribution is to show the convenience of using modern 
computer programming techniques when applied to the Dynkin approach of algebra 
representation theory. Indeed, we have been able to construct a versatile algebra- 
manipulation package, named WBasep such that: 1) WBase is a compact project, 
and the memory the Dynkin approach needs is used at best so that it works also 
on small computers; 2) WBase is easily upgradable with features. 

With regard to the point 2) in this work we will describe the new features 
of the WBase V2.0 version, a) the new table-generation routines; and b) the new 
Young-display support routines. The WBase V2 . now supports directly all Cartan 
algebras, both classical and exceptional. 

Algorithms more specialized but faster than the Dynkin one use the extended 
Young diagrams. In WBase V2 . we added the extended Young diagram display 
capability for the classical algebras in order to analize these alternative methods. 



*Supported in part by M.P.I. This work is carried out in the framework of the European Community 
Programme "Gauge Theories, Applied Supcrsymmetry and Quantum Gravity" with a financial 
contribution under contract SC1-CT92-D789. 
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2. Dynkin's approach to representation 



theoryQ 



A simple Lie algebra in the Cartan-Weyl basis is described by a set of I simultane- 
ously diagonalizable generators Hi and by the other generators E a , satisfying 

[H i ,H j ]=0 i=l,. ..,1 (1) 

d — I d — I 

[Hi, E a \ = otiE a , i = 1, a = — ; (2) 

I is the rank of the algebra, and the set of /-vectors on are the roots. It results that 
all roots can be constructed via linear combinations by a set of I roots, called simple 
roots. 

By the simple roots one then constructs the / x I Cartan matrix, which is the 
key to the classification of the Lie algebras, and it is known for all of them: the A n , 
B n , C n , D n series and the exceptional algebras G2, F4, E 6 , E 7 , E%. In the Dynkin 
approach the Cartan matrix is all we need to completely describe the algebra. This 
is at the basis of our package: the routine wstartup, given the name of the algebra, 
takes care of generating algorithmically the related Cartan matrix (called wcart in 
WBase). 

The metric Gij , which is related to the inverse of the Cartan matrix (called 
wmetr in WBase) introduces a scalar product in the space of weight vectors, which 
are Z-uples of integer numbers. Each irrep in an algebra is uniquely classified by a 
weight vector, the highest weight A whose components are all positive integers (the 
Dynkin labels) . The different states in an irrep in a given irrep are again described 
by a weight vector w; the full set of all states of a given irrep is thus described by 
a set of weight vectors, called the weight system. 

The dimension of an irrep A can be calculated with the help of the Weyl formula 
(encoded in the weyl function), 

d -( A )= n ^f 1 ® 

. (<5, a) 

pos.rootsa y ' 

where A is the highest weight determining the irrep, S — (1, ■ ■ ■ ,1), and ( , ) is the 
scalar product constructed with the metric Gy . The positive roots are the positive 
weight vectors of the weight system of the adjoint irrep. 

The weight system of an irrep is needed to reduce the products of irreps, which is 
one of the main task of our package. The weight system is obtained by subtracting 
the simple roots ai from the highest weight as described by the following recursive 
procedure (handled by wtree): 



for all i in 1 , . . . , I 

compute Uk = A — ken, k = 1, • • • ,u(i) 
add all Uk to the weight system 
restart the procedure with A := Uk 
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The computational heaviest part of the construction of the weight system is 
the computation of the degeneration of each weight vector. It is computed by the 
Freudenthal recursion formula (encoded in f reud), which needs the degeneration of 
previous levels, 



deg(u>) 



2 J2pos.roots a Efc>o dc S<> + ka) (w + ka, a) 



A + S 



w + S 



(4) 



with the initial condition deg(A) = 1. 

Once developed the machinery to generate full weight systems along with their 
degenerations, we can reduce the products of irreps, i.e. solve for the A; in the 
equation 

A a ®A 6 = A 1 ©---eA„, (5) 

each with its degeneration, where each highest weights A identify uniquely an irrep. 
The algorithm (encoded in wpdisp) is the following, wsyst: generate the weight 
systems ws(A a ) = {w an }, ws(Ab) — {w& m }, and (bprod): construct the set P = 
{uJan + Wbm}\ then, until P is empty, (whighest): find the highest weight A, in 
P, (wsyst): generate the weight system los(Aj) = {wik}, and (bremove): subtract 
from P each weight vector in los(Aj). 

3. Dynamic data structures 

The fundamental data type underlying our routines is the dynamically allocated (by 
walloc/wf ree) weight vector wvect, with a run-time chosen length wsize = I. The 
extended weight vector wplace contains also the degeneration and level informations 
of the weight vectors embedded in the weight systems. The wplace data type is 
used essentially to give a structure to the raw data kept in the blocks on which we 
base our granularity-flexible allocation scheme for the weight systems. 

The difficulty in constructing such weight systems for arbitrary irreps is due to 
the fact that the dimension of the table needed to store the weight systems is not 
known in advance. The only solution that does not waste large amounts of memory 
and does not limit our routines more than the hardware does is a multiple block 
allocation scheme. The data structure is as follows: 





a singly linked list of blocks of identical size fixed at run-time (according to the 
fragmentation required) by the number of wplace vector entries as given in bsize; 
the typedef that defines a single block leaves therefore its main structure undefined, 
as it is different for different ranks wsize and block dimension bsize: 
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typedef struct tblock { 
struct tblock *next ; 
char body[l] ;} wblock; 

This definition, as the wplace one, is used through casts that give form to unstruc- 
tured raw data as returned by the allocator balloc. Single blocks are released with 
a call to bf ree, linked blocks are released with a call to bsf ree for the first of them. 
The structured form of the blocks, when casting with wplace and defining bsize 
and wsize, is as follows: 
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4. High level routines 

The task of setting up the algebra structures and handle the lists of vectors is taken 
by the following procedures: 

wstartup (type , rank) , wcleanup (string) 
These are the open/close routines for setting up the algebra; given e.g. the algebra 
B5 (SO(ll)) one would call wstartupOB' ,5) to allocate the related vectors and 
matrices, and then wcleanupO to deallocate them. Each invocation of wstartupO 
must be followed by wcleanupO. Multiple invocations of this pair of routines are 
needed for the construction of a table of several algebras. 

wread(hw), wfdisp(hw), wydisp(hw) 
These are the input/output routines, wread read a highest weight vector label from 
the standard input, wfdisp(hw) calls the functions weyl, casimir and wheight to 
display the related irrep informations on the standard output, wydisp is the new 
WBase V2 . extended Young tableaux display function available for all non excep- 
tional algebra types. These routines need a valid wvect as handled by walloc ()/wf ree () . 

wsave (w,base , level) , wremove (w,base) 
Insert /remove the weight vector w in the list pointed by base. If level is greater 
than zero, wsave does a sorted insertion without degeneration increment, otherwise 
w is stored with its degeneration, wremove will remove a weight vector w from the 
list base only if the degeneration is just 1. 

wsyst(hw), bdisp(base), bsf ree (base) 
wsyst returns the full linked list of blocks of the weight system of highest weight hw 
along with the degenerations. It uses the recursive algorithm described in section 
2, by calling wsave to store the weight vectors, then computes the degeneration 
with the Freudenthal formula, bdisp(base) displays to the standard output all 
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entries of the weight system in base. The block list constructed by wsyst has to 
be deallocated with bsf ree. 

wpdisp (hwl , hw2 , mod) 
This function hides all the complexities of reducing products of irreps and of the 
underlying data structure, by giving to the standard output all the irreps in the 
product, from the highest to the lowest, according to the modality choosen by mod? 
The product routines are also available as iteration functions (see below). 

User interface in WBase V2.0 
In the file wmain . c it is implemented an ANSI C terminal- like interface with the 
user. A more sophisticated user interaction may be constructed taking this file as 
an example. Thanks to the new capabilities introduced in WBase V2.0, we had to 
add more options which are still one-letter options (for the details, refer directly to 
the source code). 

5. Iterators 

One of the standard conceptual device of the object-oriented technology is the iter- 
ator. The iterator hides the implementation details providing a consistent interface 
for moving through the data structure of the object considered. In WBase we in- 
troduced: 1) the iterator needed to scan the wblock list which contains the weight 
tree and 2) the iterator used to generate the wblock lists which contain the decom- 
posed irreps of a product. New in WBase V2.0 is 3) the table-generation iterator. 
Following is a short description of their use. 

Scanning through the wblock list 
It is done through the pstart/pnext iterators: 
wplace *p; 
pcurr pc; 

f or (p=pstart (base ,&pc) ; p ! =NULL ; p=pnext (&pc) ) 
< do something with p > 

remembering that: p->vect gives the weight vector, p->deg its degeneration, and 
p->level the level of the vector within the weight system. To remove the last entry 
from the list one uses base=plast (p,base), with the just removed vector returned 
in the area pointed by p. 

Getting the irreps of the decomposition of products 
The iteration functions wpstart (base) /wpnext (base ,b) that interface the con- 
struction of products return a wblock pointer to the full weight system of the 
reduced irrep. In order to display the weight tree of each irrep in the product of 
the two irreps with highest weight hwl and hw2 the fragment of code is as follows: 
wblock *base,*b; 

base=bprod (wsyst (hwl) ,wsyst(hw2)) ; 
f or (b=wpstart (base) ;b ! =NULL ; b=wpnext (b,base) ) 
bdisp (b) ; 



a See the header file wbase.h. 
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Generating irreps of increasing dimensions 
This is one of the main innovations of WBase V2 . that were only announced when 
we first introduced the package □. The table-generation procedure wsequence is 
used to produce increasing dimension irreps: 
wblock *b=NULL; 
wvect hw=walloc(); 
int dim; 

while (b=wsequence (hw , b , maxdim , fedim) ) 
wf disp(hw) ; 

wsequence provides to: 1) allocate the first block when invoked with b=NULL the 
first time, 2) allocate eventual subsequent blocks, to store the encountered highest 
weights hw, and 3) deallocate all of them when maxdim is reached. In this last 
case it returns NULL. Like the other iteration procedures, wsequence has as an 
argument a specific iterator object (the wblock pointer b) which identify each irrep 
sequencing. It is then possible to nidificate multiple iterators to produce table of 
products. 

In the next version of WBase we will probably transform the algebra initializa- 
tion/destruction functions wstartup/wcleanup in iterators in order to extend the 
table-generation capabilities of WBase V2 . to multiple algebra tables. We will also 
probably switch to the CH — h language to provide a more consistent iteration in- 
terface across the different data structures. The operator overloading capabilities 
of the CH — h language will also simplify the interface of our package by unifying 
our list storage/remove functions (by overloading the += and -= operators) and our 
input/output functions (by overloading the << and >> operators). 

6. Conclusions 

The WBase package, in our opinion, represents a useful and self-contained demon- 
stration of the convenience of new object-oriented software technology when com- 
bined with the C powerful dynamic data allocation facilities. As a by-product, we 
obtained what we think can be a useful tool for not too heavy irrep-related necessi- 
ties of physicists. The Dynkin approach of Lie algebra representation theory helped 
us to maintain a unified and elegant structure to our .package, however it must be 
noted that there are less general but faster algorithmsa based on tensor and spinor 
manipulations useful in computing products of irreps. 
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